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This paper states that as the temperature tends to the degeneracy temperature, 
there are sufficiently many particles accumulating at the point of Bose condensate 
and a gap is being formed in the spectrum of the Schrodinger equation corresponding 
to the Bose gas. In the two-dimensional case, the spectrum is purely discrete, there 
is no transition to the integral formula, and the gap in the spectrum becomes more 
visual [12j. The case of finitely many particles of the order of 10 2 is considered in 
the paper on the basis of heuristic considerations similar to those used by V. A. Fock 
to derive the Hartree-Fock equation. The paper is written in the physical language. 
The author's old results about the A-point and superfluidity are given. The author's 
paper written in 1995 is presented in the Appendix. 



In 1925, Einstein, when examining a work of Bose, discovered a new phenomenon, 
which he called the Bose condensate. A modern presentation of this discovery can be 
found in [lj. An essential point in this presentation is to define the entropy of the Bose 
gas. The definition is related to the dimension by means of the so-called "number of states" 
(cells), which is denoted by Gj in the book [lj. After this, the problem of minimizing the 
entropy is considered by using the Lagrange multipliers under two constraints, namely, 
i ! for the number of particles and for energy. The number of states Gj is determined by the 
formula which mathematicians call the "Weyl relation;" it is described in detail in [2] in 
the "semiclassical case" in the section "Several degrees of freedom." The 2-D-dimensional 
phase space is partitioned into a lattice, and the number Gj is defined by the formula 



1 (2tt/i) 



D 



The indeterminate Lagrange multipliers are expressed in terms of temperature and chem- 
5^ ! ical potential of the gas. 

Further, in [T], following Einstein, a passage to the limit is carried out as A" — > oo, 
which enables one to pass from sums to integrals. Then, in the section "Degenerate Bose 
gas," a point is distinguished which corresponds to the energy equal to zero. This very 
point is the point of Bose condensate on which excessive particles whose number exceeds 
some value A^ ^> 1 are accumulated at temperatures below the so-called degeneracy tem- 
perature Trf. The theoretical discovery of this point anticipated a number of experiments 
that confirmed this fact not only for liquid helium but also for a series of metals and even 
for hydrogen. 

From a mathematical point of view, distinguishing a point in the integral is an incorrect 
operation if this point does not form a S function. In particular, for the two-dimensional 
case, this incorrectness leads to a "theorem" formulated in various textbooks and claiming 
that there is no Bose condensate in the two-dimensional case. 



1 



In this paper, we get rid of this mathematical incorrectness and show that, both in 
the two-dimensional and in the one-dimensional case, the Bose condensate exists if the 
point introduced above is well defined. 

The main idea of the author in the proof of the occurrence of the Bose condensate in 
the D- dimensional case is in a concordance between the chemical potential \i — > and 
the number of particles iV — > oo when passing to the limit. 

The phenomenon associated with the point of condensation holds only if the limit as 
\i — > depends on N — >■ oo. 

If we accept Einstein's remarkable discovery for the three-dimensional case and justify 
it in a mathematically correct way, then the Bose condensate in the two-dimensional 
case is equally correct mathematically. We dwell on the two-dimensional case below in 
particular detail. 

Thus, we consider the case in which N ^> 1, but n is not equal to infinity. In the 
section "Ideal gas in the case of parastatistics" of the textbook by Kvasnikov [3J, there is a 
problem (whose number in the book is (33)) which corresponds to the final parastatistics 

1 k+1 Nj 

n > ~ exp{ £^} _ 1 " exp{(A; + l)^}-l' " Gj' U 

In our case, we have k = Nd, and the point of condensate is Eq = 0. 

By ([I]), it is clear that Gj is associated with the D-dimensional Lebesgue measure 
and, in the limit with respect to the coordinates Aqj, gives the volume V in the space of 
dimension 3 and the area Q in the space of dimension 2. The passage with respect to the 
momenta Apj is also valid as iV — > oo and fi > 5 > 0, where 5 is arbitrarily small. 

Expanding ([2]) at the point Eq = in the small parameter 

x = (fiN d )/T d , 

where N d stands for the number of particles corresponding to the degeneration and T d for 
the degeneracy temperature, and writing 

i = -n/T d , 

we obtain Go = 1, see ( TT2"j) below, 

N d + 1 \ _eW*-l- (N d + l)(e« - 1) 

no 



exp{^} - 1 exp{(iV d + !)=£} - 1 j (e* - l)(e 2 ^ 
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x 2 - 0.19Lr d - . . . . (3) 



For example, if x — » 0, then n$ = N d /2, and hence the number no in the condensate 
at T = T d does not exceed Nd/2. If x = 1.57, then n$ m N d /10. Certainly, this affects the 
degeneracy temperature, because this temperature can be expressed only in terms of the 
number of particles above the condensate, iVd, rather than in terms of the total number 
of particles Nd (which is equal to the sum of N d and of the number of particles in the 
condensate) . 

According to the concept of Einstein, at T = Td the condensate contains o(N d ) parti- 
cles. However, even this accumulation gives a 5 function, albeit with a small coefficient, 
for example, N d /hiN d . 



To reconcile the notion of Bose statistics which is given in |lj with symmetric solutions 
of the A-particle Schrodinger equation, i.e., of the direct sum of iV noninteracting Hamil- 
tonians corresponding to the Schrodinger equation, and the symmetric solutions of their 
spectrum, it is more appropriate to assign to the cells the multiplicities of the spectrum 
of the Schrodinger equation in the way described in j3]. 

Consider the nonrelativistic case in which the Hamiltonian H is equal to 

p 2 /(2m), 

where p stands for the momentum. 

The comparison of Gi with the multiplicities of the spectrum of the Schrodinger equa- 
tion gives a correspondence between the eigenf unctions of the iV-partial Schrodinger equa- 
tion that are symmetric with respect to the permutations of particles and the combina- 
torial calculations of the Bose statistics that are presented in pQ. 

A single-particle -^-function satisfies the free Schrodinger equation with the Dirichlet 
conditions on the vessel walls. According to the classical Courant formula, 



X j ~ — ( -77 ) T as j ->■ oo, (4) 



2hW n D / 2 T(D/ 2 + l)\ 2/D 
m \ V 

where D stands for the dimension of the space, because the spectral density has the 
asymptotic behavior 

y m D/2^D/2 

p(A) = rWTWw (1 + 0(1)) as A ^°°' (5) 

The asymptotics (pjj is a natural generalization of this formula. 

Using this very correspondence, we establish a relationship between the Bose-Einstein 
combinatorics [T], the definition of the A-particle Schrodinger equation, and the multi- 
plicity of the spectrum of the single-particle Schrodinger equation. 

The spectrum of the single-particle Schrodinger equation, provided that the interaction 
potential is not taken into account, coincides, up to a factor, with the spectrum of the 
Laplace operator. Consider its spectrum for the closed interval, for the square, and for 
the D-dimensional cube with zero boundary conditions. This spectrum obviously consists 
of the sum of one- dimensional spectra. 

On the line we mark the points i = 0, 1, 2, . . . and on the coordinate axes x, y of the 
plane we mark the points with x = i = 0,1,2, .. . and y = j = 0,1,2, ... . To this set of 
points we assign the points on the line that are positive integers, I = 1,2, ... . 

To every point we assign a pair of points, i and j, by the rule i + j — I. The number 
of these points is nj = I + 1. This is the two-dimensional case. 

Consider the 3-dimensional case. On the axis z we set k — 0, 1, 2, ... , i.e., let 

i + j + k = I 
In this case, the number of points n\ is equal to 

(Z + l)(/ + 2) 

m = . 

2 

It can readily be seen, for the D-dimensional case, that the sequence of multiplicities 
for the number of variants 



D 

I = 

k=l 
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where rrik are arbitrary positive integers, is of the form 

*M = (i-twP-1)! ' forjD = 2 ' *( 2 ) =< . ( 6 ) 

OO OO 

J2n* = N, eY,Qi{D)Ni = E. (7) 

i=l i=l 

The following problem in number theory corresponds to the three-dimensional case 
£> = 3(cf. [lj): 

f> = iV, ^^p.N t = E, - = M. (8) 

i=l i=l 

Write M = E d /ei, where E\ stands for the coefficient in formula (0} for j = 1. Let us 
find E d , 

f°° l -fde 

Ed= A ? (9) 
./o e~< d — 1 



where 

bl 2 <#>i . . . dpp dV D 
2m (2irh) 

Whence we obtain the coefficient a in the formula 



de = — 7KZTT5 • ( 10 ) 



E d = aT 2 /\(l + D/2h(l + D/2), (11) 

where 7 = D/2 — 1. 

To begin the summation in (J7|) at the zero index (beginning with the zero energy), it 
is necessary to rewrite the sums in the form 

Y^N l = N, e y ST(q i (D)-l)N i = E-eN. (12) 

i=0 i=0 

The relationship between the degeneracy temperature and the number Nd of particles 
above the condensate for fi > 5 > (where 5 is arbitrarily small) can be found for D > 2 
in the standard way. 

Thus, we have established a relationship between Gi in formula ([T]) (which is combina- 
torially statistical) and the multiplicity of the spectrum for the single-particle Schrodinger 
equation, i.e., between the statistical [1] and quantum- mechanical definitions of Bose par- 
ticles. 

For D = 2, the general problem reduces to a number theory problem. 
Consider the two-dimensional case in more detail. There is an Erdos' theorem for a 
system of two Diophantine equations, 

00 00 

j2n 1 = n, Y, tN * = M - ( 13 ) 

1=1 8=1 

The maximum number of solutions of this system is achieved if the following relation is 

satisfied: 

N d = c^M 1 / 2 log M d + aM x J 2 + o(M d 1/2 ), c = vr y^2/3, (14) 
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and if the coefficient a is defined by the formula 



c/2 = e 



-co/2 



The decomposition of into one summand gives only one version. The decompo- 
sition Md into Md summands also provides only one version (namely, the sum of ones). 
Therefore, somewhere in the interval must be at least one maximum of the variants. Erdos 
had evaluated it (ITS]) (see [7]). 

If the number N increases and M is preserved in the problem (Tl3|) . then the number 
of solutions decreases. If the sums ( Tl3l) are counted from zero rather than from one, i.e., 
if we set 

oo oo 

Y^m = (M-N), ^iV, = iV, (15) 

i=0 i=0 

then the number of solutions does not decrease and remains constant. 

I'll try to explain this effect. The Erdos-Lehner problem [5j is to decompose Md into 
N < Nd summands. Let us expand the number 5 into two summands. We obtain 3 + 2 = 
4+1. The total number is 2 versions (this problem is known as "partitio numerorum"). If 
we include to the possible summands, we obtain three versions: 5 + = 3 + 2 = 4 + 1. 
Thus, the inclusion of zero makes it possible to say that we expand a number into k <n 
(positive integer) summands. Indeed, the expansion of the number 5 into three summands 
includes all the previous versions, namely, 5 + + 0, 3 + 2 + 0, and 4 + 1 + 0, and adds 
new versions, which do not include zero. 

In this case, the maximum number of versions for the decomposition of the number 
5 into iV summands (there are two versions) is achieved at N = 2 and N = 3 (the two 
values for the maximum number of versions for N above the condensate). 

In this case, the maximum does not change drastically [5]; however, the number of 
versions is not changed, namely, the zeros, i.e., the Bose condensate, make it possible 
that the maximum remains constant, and the entropy never decreases; after reaching the 
maximum, it becomes constant. This remarkable property of the entropy enables us to 
construct an unrestricted probability theory in the general case [8]. 

Let us turn to a physical definition. 

Note first that, without changing the accuracy of the quantity whose logarithm is 
evaluated, we can replace log Md by 

(l/2)log(iVQ). 

Then 

2N d /Q fNd\ 



c- 1 \og(N d /Q) + a \Q 

In our case, Nd/Q corresponds to the number of particles above the condensate. 

According to formula ( II lj) . in the two-dimensional case we must set 7 = and find 
the coefficient a. Then formula ( 1T61) gives us a relationship between Nd and Td due to the 
fact that the number of particles in the condensate is o(Nd). 

For details concerning the Bose condensate in the one-dimensional case, see j 12 j . 
and [T3] . 

In fact, we have proved that there is a gap between E\ > and Eq = for \x = 0. 

Remark 1. Transition to quantum perfect liquid occurs for T = 0.84Td- For 0.84 < 
T < Td, there is a two-phase region up to the point of equality of chemical potentials. 
But this liquid is not superfluid. As follows from the above and as was noted earlier 



by N. N. Bogoliubov, some " weak imperfectness" is required for the superfluidity. In 
this case, the liquid spinodal in the region of negative energy terminates at the point 
T = 0.39Trf. Below this temperature, the liquid becomes a dispersive system, i.e., "foam", 
cf. larXiv:1111.6l"06 v2. 

The liquid spinodal 7 = 7(T) for 7 < is determined by the following conditions: 

iV = T r 7 0+1 C(7o + l), 7o = l/2, T r = T/T d , 

on the interval Z = Zj, where Z = PV/NT, in the plane Z, P; for the isotherm-isochor 

77 0+1 (C(7o + l)) = C(l + l7(T r )|)T r , 

we obtain the dependence of 7 on the temperature for fl = (for the spinodal) in the 
region of negative pressures. In particular, for T r — 1 (i.e., for T = T^), 7 = —70 
determines the initial reference point for the chemical potential of the liquid phase jx = 
by the condition 

p _ Li 2 _ 70 (e-^) _ i 
C(2 + 7o) 

Thus, we find the initial reference point jx from the condition on T r : 

Li 2+7o (l) = Li 2 „ 70 (e-^), 
where Li is the polylogarithm. We find 11 = jl as a function of Td from the condition 

77 Li 2+70 (e"^) = T"l^)l Li 2H7(Tr) |(e-^^/^). 

This allows us to obtain the point \i on the gas branch, which determined the phase 
transition of the perfect Bose gas into the perfect Bose liquid. The two-phase region 
is possible only on a small interval T r < 0.84. For lesser temperature T r , the chemical 
potential/i = —00, i.e., the gas completely turns into liquid. Only in the case where the 
interaction is taken into account (see below), there is a sufficiently good coincidence with 
experimental data. 

Remark 2. The author studied the relationship between the economy during a cri- 
sis and the Bose condensate, which corresponds to the bankruptcy |9j. Continuing the 
correspondence principle proposed by Irving Fisher, an economist and a disciple of Gibbs 
(this principle is the "fundamental law of economics"), where the amount of money M 
corresponds to the number of particles N, the author suggested to compare the chem- 
ical potential to the negative value of the nominal interest rate, which corresponds to 
Friedman's rule. 

The issue of money accompanied the fall of the nominal interest to 0.5% following this 
dependence in which the small parameter jrNd became equal to ^75, where D stands for 
the "number of degrees of freedom", which can be fractional (in our case, this number is 
the dimension) |8j. 

In the paper [10J by E. M. Apfel'baum and V. S. Vorob'ev, taking into account the de 
Boer parameter, the values Z c for helium were calculated experimentally, which, according 
to the author's rule 

7 C(f + 1) 
c " C(f) ' 



enables one to determine the number D of degrees of freedom and to experimentally 
verify whether or not there is an empirical relation of this kind between \x and Td in 
thermodynamics. It was assumed that the Bose gas is not perfect (i.e., the Schrodinger 
equation with a potential is considered) and the value of x = (fiN^/T^ reflects the 
interaction between the particles, just as the Zeno line reflects the interaction between 
particles in classical thermodynamics 

Let us present heuristic considerations concerning the passage through the point Td 
which were presented by the author in [T3]. The author proved and used the Bose statistics 
in the case of a fractional number of degrees of freedom for classical thermodynamics, 
where to a value of Td there corresponds the critical temperature. The author has shown 
that these values coincide. For an example describing the creation of a dimer, it is shown 
that, for T = T c , one degree of freedom becomes " frozen", and we obtain two degrees 
of freedom rather than three. For a dimer with T > T c , if the oscillational degrees of 
freedom are taken into account, then the number of degrees of freedom becomes equal 
to 6. Two degrees of freedom are obtained under the assumption that the oscillational 
degrees of freedom of the dimer are also "frozen" at T = T c . If we suppose this heuristic 
supposition for the quantum case, then, for T < T^, both dimers with two degrees of 
freedom and dimers with six degrees of freedom are created. This corresponds to the 
two-liquid Thiess-Landau model. In this case, the dimers with two degrees of freedom 
give the A-point and the dimers with six degrees of freedom give superfluidity. Indeed, in 
the two-dimensional case we have 



and we obtain a logarithmic divergence at the point £ = for /i — > 0. 

Thus, if we consider an iV-particle Schrodinger equation whose eigenfunctions are 
symmetric under the permutations of the particles, then the parastatistic correction leads 
to the fact that N/2 particles are in the condensate for T = Td and N = N^. For N > Nd, 
all the extra particles pass into the condensate state, which determines the dependence of 
the temperature T on N, and hence the dependence of iV on the temperature for T < Td 
as well. 

The case in which N is not so large as it is in statistical physics, i.e., the so-called 
mesoscopic state (see |14|). can also be of interest for us. In this case, let us use Fock's 
idea for the Hartree equations, which lead to the Hartree-Fock equations. *** Namely, we 
consider the single-particle equation of the mean field (a self-consistent field) and apply (to 
the resulting "dressed" potential) the procedure of transition to the iV-partial Schrodinger 
equation with a dressed potential, just as we proceeded above for the operator j-A. 
Here we can consider two ways of investigation. The first way is the way used by Fock 
and which leads in the semiclassical limit to the Thomas-Fermi equations for the dressed 
potential. Another way is to consider the Hartree temperature equations (see [15J) and 
to obtain the Thomas-Fermi temperature equations in the classical limit. 

Since the quantity T d is small, it is easier to use the first way and to find the "dressed" 
potential. 

Let V(q — q') be a pairwise interaction potential such that J \ V(r)\dr < oo. The 
dressed potential W(q) is given by the formula 
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where U(q) stands for the external potential and ip(q') for an eigenfunction of the 
Schrodinger equation which depends on the "dressed" potential and is thus an equation 
with a "unitary" nonlinearity. The expansion of the equation in powers of h can be found 
by the method of complex germ up to 0(h k ), where k is an an arbitrarily large numbed 
(see [23], where system (63) defines a complex germ; see also [24j-[30j). The superfluidity 
in nanotubes was confirmed experimentally. 

The author thanks Professors G. I. Arkhipov, V. S. Vorob'ev, and V. N. Chubarikov 
for permanent discussions. 
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Quasi-Particles Associated with Lagrangian Manifolds 
Corresponding to Semiclassical Self- Consistent Fields. Hid 
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In the preceding part of this paper, we presented Eqs. (25) for quasi-particles associated 
with an n-dimensional Lagrangian manifold and Eq. (29) for quasi-particles corresponding 
to a (2n — l)-dimensional manifold. These equations were written out only in the x-chart, 
and the quantum corrections were given without proof. In this part we essentially use 
the canonical operator ideology to obtain Eq. (25) with corrections in the x-chart as well 
as in any other chart of the canonical atlas [1]. To derive the correction in Eq. (29), a 
"modified" 5-function must be used, and this will be done in the next part of the paper. 

To obtain the result in an arbitrary canonical chart, one should pass on to the p- 
representation with respect to some of the coordinates in the Hartree equation. This is 
actually equivalent [2] to considering the Hartree-type equation 

#o( A x, - -4 ih-^-j + J dyip*(y)H x ^A- x, - -4 ih-^-; -4 y, - -4 ih^-jip(y) ip{x) 
= Qip{x), (50) 

where x,y G lR n , ip G L 2 (M. n ) is a complex-valued function, h > 0, Q G M, and the indices 
1 and 2 specify the ordering of the operators x and —ihd/dx. The function Hi satisfies 
the condition Hi(x,p x ;y,p y ) = Hi(y,p y ; x,p x ). Equation (50) generalizes the ordinary 
Hartree equation (Eq. (1) in [4], where iV = 1). The study of Eq. (50) is important, for 
example, if one makes an attempt to find a solution to the Hartree equation (1) in the 
momentum representation, 



ip(x) = / ip(p)e 



,(i/h)px dp 



Let us also discuss the variational system associated with Eq. (50), which can be 
obtained as follows. Along with Eq. (50), let us write out the conjugate equation and 
consider the variations of both equations assuming that the variations Sip = F and Sip* = 
G are independent. 

The variational system has the form 

ffo(Ai,-A ih ^) ~ Q + J d y^*(y) H i( ^ ih-^;Ay,- 4 ih^jip(y) F(x 

2 The work was supported by the ISF under grant No. MFO000. 



in 



+ J dy(G(y)H 1 ^ x,- 4 ih ^A V,- 4 ih^jip(y) 

+ r{y)H, ( 4 x, - 4 ih— ; 4 y, - 4 i/i— )F(y))^(x) = -/3F(x), 

#o( 4:r,4 i/ij^) -Q + J dy^H^Ax^ih^Ay^ih-^j^iy) G(x) 

dy^H, ( 4 x, 4 i/i— ; 4 y, A ih—)r{y) 

+ ( 4 4 4 y, 4 ih—)G(yj)r(x) = PG(x). (51) 

Equations (50) and (51) play an important role in the problem of constructing asymptotic 
solutions to the iV-particle Schrodinger equation as N — > oo [5-7]. 

For example, the spectrum of system (51) (possible values of j3) corresponds to the 
spectrum of quasi-particles. Namely, the difference between the energy of an excited state 
and the ground state energy is given by the expression 



where the numbers G Z + , k = l,oo, which are equal to zero starting from some k, 
define the eigenfunction and the eigenvalue of the excited state, and (3k G K. are the 
eigenvalues of system (51). 

In this paper we are interested in asymptotic solutions to Eqs. (50) and (51) as the 
"inner" h tends to zero. 

Asymptotic solutions to Eq. (50) are given [8] by the canonical operator on a La- 
grangian manifold A™ = {x = X(a),p = P{ct)} invariant with respect to the Hamiltonian 
system 

^ = dH(x,p) . = dH(x,p) , 52 
dp dx 



where 



H(x,p) = H (x,p) + J d/j, a H 1 (x,p;X(a),P(a)), 



a G A™, and d\i a is an invariant measure on A n . The Lagrangian manifold lies on the 
surface H(x,p) = Q. If a chart A is projected diffeomorphically in the rr-plane, then 
the canonical operator acts as the multiplication by exp{(i/h)S(x)}/^/j, where S(x) = 
J pdx on A™ and J = DxjD\i a . We are interested in finding asymptotic solutions to 
Eqs. (51). Without loss of generality, we can consider only the case of a>chart. Indeed, to 
obtain similar expressions in the p-chart, one must consider the Fourier transformation of 
Eqs. (50) and (51) and apply the same technique, since the form of the equations remains 
unchanged. 

Let us seek the asymptotic solutions to Eqs. (51) in the x-chart in the form 

F(x) = f(x)^(x), G(x) = g(x)P(x), (53) 

where the functions / and g, in contrast to ip and ip*, have a limit as h — > 0. One can 
consider a more general case, by allowing / and g to be functions of x and —ihd/dx, but 
in the leading term as h — > we have 

ox ox 



1 1 



and so we arrive at functions / and g that depend only on x. 

The second equation in system (51) can be rewritten in the form 

+ j dy\y>(y)Ay)H 1 (x,ih^;y,ih^j'ijj*(y) 



//| (•'••'/' 77- ://•/'/' 37-): 

ox oy> 



(yj\r{y)} = Pg(x)r(x), 



(54) 



where [A; B\ — AB — BA and 

c(x) = f(x)+g(x). (55) 

Equation (50) is used in the derivation of Eq. (54). We observe that all terms containing 
the function g on the left-hand side in Eq. (50) are 0(h), since the commutator of two 
operators depending on x and —ihd/dx is equal in the classical limit to (—ih) times the 
Poisson bracket of the corresponding classical quantities. 

Thus, the function c, as well as the eigenvalue /3, is assumed to be 0(h). Let us rescale 
these quantities as follows: 

c( x ) = hc(x), /3 = h/3. (56) 

Now we can derive the equation for g,c, and (3 in the leading term in h and the first 
correction to it from Eq. (54), making use of the following relations: 



} [ A ( 4 *■ 4 ih iL) ; = J ih Wa ( 4 4 ih m) 



d \ w j ., OA / 1 2 .,d\d£ 

a \ Ox) Ox a 

h 2 d 2 A ( 1 2 ., 9\ <9 2 £ 
— )> rr, — > ih — I 



a, o=l 



(9x7 dx a dxb 



where p a = ihd/dx a , A(x,p) is a function M 2n ->• C, £ : R n ->• R. 

ii) ip(x) = x( x : h)e^ h ^ x \ where x — l/V~J in the leading term in h; 



iii) 



ox V Ox„ ox J 



(57) 



iv) 



/ , d dS\ „(dS\ 9 ih ^ d 2 B d 2 S 



-dS\ 



+ 



+ 



(ihf 



E 



d 2 B d 2 



f-f dp a dx a 2 4^ <9p a <9p <9x a &r 

a=l a, 6=1 

i/i) 2 ^ <9 3 5 <9 2 ,S 9 



2 ^dp a dp b dx a dx b 2 j-^ dp a dp b dp c dx a dx b dx c 

a, o=l a,o,c=l 



E 



5 ^ dp a dp b dpc dx a dx b dx c 

a,b,c=l 

d*B 2 S d 2 S 



8 t j , dp a dp b dp c dp d dx a dx b dx c dx d 

a,b,c,a=l 

where all derivatives of B are evaluated at the point p = dS/dx. 



+ 0(h 3 ), 



(58) 
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These relations can easily be obtained for monomial functions A and B. An application 
of formulas i)-iv) yields the equation 

+ t j ^ _ + _ yj (A - (a)) __ (A>)) 

— X Cli^O — _1_ 



2 ^ <9x a 

a,6,c=l 



+ | / ^c(X(/3)) £ [ 



dPa^Pb^Pc dxbdx c 
d 2 H l d 2 S 



tdpldpt dx a dx b 



(*(«)) + 



2 ^ dx a dx b 

a, 6=1 

d 2 ^ d 2 £ 

dp y a dp\ dx a dx b 



dPadpl 



W))] 



ih 
~2 



^ L dp* &r a 



<9pa <9x a 



+ 



a=l 6=1 



<9 2 #i 9 In J 



d 2 H l din J 



6,c=l 

2 



<9 2 S 



dpldpfdpl dx h dx, 



{X(a)) + 



dpldpl dx b dpldpl dx b 

d 3 ^ d 2 S 



(X(a)j) 



dpldp\dp v c dy b dy t 



W))} 



/<9 

a, 6=1 



8 2 H X d 2 g 



dpi dx a dx b 



W)) = 0; 



(59) 



in this formula the arguments 

x = X(a), p x = P(a), y = Xtf), p y = P((3) (60) 

of the function H 1 and of its derivatives, as well as the arguments x = X(a), p x = P(a) 
of the function H, are omitted. 

Let us now find another equation relating g to c. To this end, let us multiply the first 
equation in system (51) by ip*(x) and the second equation by ip(x). Let us subtract the 
first product from the second. We obtain 

(x,ih^);g(x) 



/3ip*(x)ip(x)c(x) = ip(x) 



Hi 



4>*(x) 



+ 
+ 



d 



ip*(x) H\x,—ih—j;g(x) ip(x) — ip*(x) H yx, —ih—j; c(x) ip(x) 
ijj(x) J dyip(y) H 1 (^x,ih^]y,ih-^-y,g(y) ip* (y)ip* (x) 

dy*P*(y) #(y);#i(a;, -i/i— -i/i— ) ij;(y)ip(x) 
+ ip(x) j dyip(y)c(y)H 1 (x,ih^;y : ih-^j'ip*(y)'ilj*(x) 

-ip*(x) J dy^*(y)H 1 (^x,-ih^;y,-ih-^c(y)^(y)^(x). (61) 
Let us use Eqs. (57)-(59). We find the following equation for g and c modulo 0{h 2 ): 



•ES^W) - - t SI (Ma)) 



dp x a dx a 



a 6=1 ®Pa®Pb dx a dx b 
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If 



Ed 2 H dg , v , ,,<91nJ\ . 



, ').r„ \dx a dx b dpldpfdpZ dy a dy b dpldp y b dp y c 

+ J ** £ ^(fl) ( ^ («)) - ^ + ^r(- Y (»'^» 



4^(0) g (_(*(„))_ + 



+ 2 ^ a^? {a)) dp x a dptdp* c dx b dx} {a)) 

a,b,c=l 

dc d 2 H d\nJ h >A <9 2 iJ <9 2 c 

"2^^? [a)) d P x a dpi dx b [ {a)) + 2 ^ dp x a dpt dx a dx b { {a)) 

*-(*(ffl)^W»))&=0. (62) 



<9a; a <9rr b dpadpl 



H (x,p x ) = p 2 x /2 + U(x), H x {x,p x ;y,p y ) = V(x,y), 
then Eqs. (59) and (62) become much simpler and acquire the form 

(iVSV - + y V^(x, X(a'))c(X(a')) <W + |(- A <? + V In JVg) = 0, 

(iVSV -/3)c- A# + Vln JV#- -(-Ac + V In JVc) = 0. (63) 

From Eqs. (63) one can approximately find the functions F and G, which are important 
for constructing approximate wave functions in the iV-particle problem as N — > oo [5]. 

Let us now relate the obtained results to the solution to variational equation for the 
Vlasov equation, obtained in the preceding part of this paper [3] . 

Let p be the projection on the function ip. Its kernel is p(x,y) = ip(x)ip*(y), and its 
symbol is p(x,p) = ip{x)ip*{p)e^/ h ^ px . The operator p satisfies the Wigner equation, which 
reduces to the Vlasov equation as h — > 0. The operator a with the kernel F(x)ip*(y) + 
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ip(x)G(y) is equal to 

v = fp + pg (64) 

and satisfies the variational equation to the Wigner equation ,^which is reduced to the 
variational equation for the Vlasov equation (20). In Eq. (64) / and g are the operators 
of multiplication by the functions / and g. We see that in the semiclassical approximation 
the symbol of a is 0(h), since 

° = [p;g} + hep 

and ^ 

a(x, P )c,H(-i^^-(x, P )^- + cp). 

Since p is the (^-function in the semiclassical approximation [3], the function a is actually 
the sum of the (^-function and its derivative. Equations (63) are consistent with Eqs. (24), 
obtained in [3] for the coefficients of 5 and 5'. Thus, the approach suggested in this part 
of the article allows us to find an asymptotic formula for a as well. 

The author is deeply grateful to O. Yu. Shvedov, whose assistance in carrying out all 
computations was invaluable. 

References 

[1] Maslov. V. P. Theorie des Perturbations et Methodes Asymptotiques (Dunod, Paris, 
1972). 

[2] Maslov. V. P. "Equations of self-consistent field", in Sovremennye problemy matem- 
atiki (1978), Vol. 11, pp. 153-234. 

[3] Maslov. V. P. "Quasi-particles associated with Lagrangian manifolds corresponding 
to classical self-consistent fields, II," Russian J. of Math. Phys. 3 (1), 123-132 (1995). 

[4] Maslov. V. P. "Quasi-particles associated with Lagrangian manifolds correspondijng 
to classical self-consistent fields, I," Russian J. of Math. Phys. 2 (4), 528-534 (1995). 

[5] Maslov. V. P. and Shvedov, O. Yu. "Quantization in the neighborhood of classical 
solution in the iV-particle problem and superfluidity," Theoret. and Math. Phys. 98 (2), 
181-196 (1994). 

[6] Maslov. V. P. and Shvedov, O. Yu. "Complex WKB- method in the Fock space," 
Dokl. Ross. Akad. Nauk 340 (1), 42-47 (1995). 

[7] Maslov. V. P. and Shvedov, O. Yu. "Large deviations in the many-body problem," 
Matem. Zametki 57 (1), 133-137 (1995). 

[8] Maslov. V. P. Complex Markov Chains and Feynman Path Integral (Nauka, 
Moscow, 1976). 



15 



